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ABSTRACT 

The recent discovery by Cowie et al. (1995) and Tytler et al. (1995) of metals in the Lya clouds shows that 
the intergalactic medium (IGM) at high redshift is contaminated by the products of stars and suggests that 
ionizing photons from massive star formation may be a significant contributor to the UV background radi- 
ation at early epochs. We assess the validity of the stellar photoionization hypothesis. Based on recent compu- 
tations of metal yields and O-star Lyman continuum (Lye) fluxes, we find that 0.2% of the rest-mass energy of 
the metals produced is radiated as Lyc. By modeling the transfer of ionizing radiation through the IGM and 
the rate of chemical enrichment, we demonstrate that the background intensity of photons at 1 ryd that 
accompanies the production of metals in the Lya forest clouds may be significant, approaching 0.5 x 10 21 
ergs cm“ 2 s" 1 Hz' 1 sr _l at z % 3 if the Lyc escape fraction is >0.25. Together with quasars, massive stars 
could then, in principle, provide the hydrogen and helium Lyc photons required to ionize the universe at high 
redshifts. We propose that observations of the He n Gunn- Peterson effect and of the metal ionization states of 
the Lya forest and Lyman-limit absorbers should show the signature of a stellar spectrum. We also note that 
the stellar photoionization model fails if a large fraction of the UV radiation emitted from stars cannot escape 
into the IGM, as suggested by the recent Hopkins Ultraviolet Telescope observations by Leitherer et al. (1995) 
of low-redshift starburst galaxies, or if most of the metals observed at 2 % 3 were produced at much earlier 
epochs. 

Subject headings : cosmology: observations — diffuse radiation — intergalactic medium — 
quasars: absorption lines 


1. INTRODUCTION 

It is widely believed that the background ultraviolet flux 
with energy greater than 1 ryd arises from quasars (Sargent et 
al. 1980; Bergeron & Stasinska 1986; Madau 1992), possibly 
augmented by other sources such as hot, massive stars in star- 
bursting young galaxies (Bechtold et al. 1987; Steidel & 
Sargent 1989) or more exotic mechanisms such as decaying 
neutrinos (Sciama 1993). This radiation is responsible for 
maintaining the ionization state of the diffuse intergalactic 
medium (IGM), the Lya forest clouds, and the metal-rich 
absorption systems. Using a model-independent argument 
from metal production, Songaila, Cowie, & Lilly (1990) sug- 
gested that the ionization state of intergalactic material could 
be attributed to star-forming galaxies if at least a few percent of 
galaxy formation occurred at redshifts z = 3-5. This possibility 
was studied in more detail by Miralda-Escude & Ostriker 
(1990). 

To date, a significant population of primeval, star-forming 
galaxies at high-z has not been detected. However, recent mea- 
surements (Tytler et al. 1995; Cowie et al. 1995; Fan 1994) 
show that a large fraction of the numerous, Lya forest clouds 
are associated with measurable C iv absorption lines. Because 
the Lya clouds have large filling factors and may contain a 
significant fraction of the baryons in the universe, the new data 
are the first evidence for a widespread distribution of metals in 
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the IGM. If QSO absorption systems trace the bulk of star 
formation occurring in galaxies at high redshifts, and if the 
initial mass function (IMF) extends to massive O-type stars, 
the stellar contribution to the local ionizing background could 
be significant. 

In this paper, we examine the contribution of massive stars 
to the ionizing background, using a new calibration of the 
Lyman continuum (Lyc) emission from O stars, accurate radi- 
ative transfer of this radiation through the IGM, and a param- 
eterization of the metal formation rate in QSO absorption-line 
systems. In § 2 we summarize the baryon content and metal- 
licity of gas in Lya forest clouds, Lyman-limit systems, damped 
Lya systems, and star-forming galaxies. In § 3 we estimate the 
luminosity of ionizing photons that accompany the metal pro- 
duction and compute the spectrum of ionizing radiation. In § 4 
we discuss the implications for the ionizing background and 
chemical evolution of the universe. 

2. HEAVY-ELEMENT ABUNDANCES AT HIGH REDSHIFT 
2.1. Lya Forest Clouds 

The term “ Lya forest ” is used to denote the numerous 
population of H 1 absorption systems observed in the spectrum 
of quasars (Weymann 1993; Bechtold 1994). With H 1 column 
densities ranging from 10 17 cm 2 down to 10 13 cm 2 and 
probably even lower, these systems are usually interpreted as 
optically thin clouds associated with the era of baryonic infall 
and galaxy formation. The clouds are observed to evolve 
rapidly with redshift between 1.7 < z < 4 (e.g., Murdoch et al. 


551 



552 


MADAU & SHULL 


Vol. 457 


1986) and typically do not show metal lines in - 1 A resolution 
spectra (Sargent et al. 1980; Tytler & Fan 1994; but see Lu 
1991). However, recent spectra at high S/N and high resolution 
obtained with the HIRES spectrograph on the Keck telescope 
(Tytler et al. 1995; Cowie et al. 1995) have shown that 50%- 
60% of the Lya forest clouds with log N H , > 14.5 have under- 
gone some chemical enrichment, as evidenced by weak, but 
measurable C iv absorption lines. The typical inferred metal- 
licities, Z Lyat , range from 0.003 to 0.01 of solar values, subject to 
uncertainties of photoionization models (e.g., Bergeron & Sta- 
sinska 1986; Steidel 1990; Madau 1991; Reimers & Vogel 
1993) and possible contributions from local sources (Giroux, 
Sutherland, & Shull 1994; Viegas & Gruenwald 1991 ; Gruen- 
wald & Viegas 1993; Petitjean, Rauch, & Carswell 1994; 
Viegas & Friaca 1995). Clearly, these metals were produced in 
stars that formed in denser gas; the metal-enriched gas was 
then expelled from the regions of star formation into the IGM. 
Although carbon need not by formed by massive stars, we 
assume that the IMF extends as a power law to O-type stars of 
mass 16-85 M 0 . 

The available data do not tightly constrain the mass density 
parameter, , of the Lya clouds. Recent spectroscopic 
observations of close quasar pairs (Bechtold et al. 1994; 
Dinshaw et al. 1994) reveal that the forest clouds have charac- 
teristic sizes much larger than previously thought, with trans- 
verse radii >200 h$Q kpc, where h 50 is the Hubble constant in 
units of 50 km s“ 1 Mpc - l . Although the modeling is typically 
done for spheres, one can generalize by considering either 
spheres of radius R or clouds with large aspect ratios, a = 
R/l = 10a lo , of transverse size R to thickness l. Assuming 
photoionization equilibrium with the diffuse UV field and 
neglecting those unlikely cases in which the line of sight enters 
the cloud through its edge, we find that the total hydrogen gas 
density of a cloud of size R = R 2 oo x 200 kpc and observed 
H i column density N H , = 10 l4 /V 14 cm - 2 is 

r H1 N Hl cos 9 ~ V /2 

o + 2 X )a A (D/J 

~ (5.6 x 1 0 “ 5 )R 2 oo 2 N } x i\ ! (a cos 0) 1 - 2 cm" 3 , (1) 

where r H } is the hydrogen photoionization rate, 6 is the angle 
between the symmetry axis of the cloud and the line of sight, 
X = «hc/"h * 1/12, J 912 = j- 2l (z) X 10 21 ergs cm' 2 s' 1 
Hz~ 1 sr 1 is the metagalactic flux at the Lyman edge, a A = 4.1 
x 10 13 T 4 0 73 cm 3 s~ l is the (case A) recombination rate 
coefficient to all levels of hydrogen, and T = T 4 x 10 4 K is the 
cloud temperature. For an ionizing radiation field with con- 
stant energy spectral index a, r H ,«(1.2 x 10 11 s~ l )J _ 2i / 
(a + 3). We adopt T 4 = 4 and a = 0.5 for the background radi- 
ation. Note that setting < a cos 6 > = 1 recovers the usual limit 
of spherical clouds. 

Press & Rybicki (1993) parameterized the observed fre- 
quency of absorbers per unit redshift by dN/dz ~ 4.2(1 + z) 2 46 
for lines of rest equivalent width greater than 0.32 A, corre- 
sponding to columns N H , > 1.5 x 10 14 cm -2 (b — 30 km s~ 1 ). 
In an Einstein-de Sitter cosmology (q 0 = 0.5), the volume- 
filling factor of this subpopulation is /% 4,2(7/ 0 R/ac cos 9) 
(1 + z) 4 - 96 , assuming that all clouds have similar sizes. Denot- 
ing with n Hl crit Ih e closure hydrogen density, we find for 
^Lya =/«H/«H.cfit Z = 3, 

n Lya *0.02 hJ J - 2 ±°° )' 12 . (2) 

\tf 10 cos 9 ) 


We will show below that massive stars associated with QSO 
absorption-line systems may generate an ionizing flux J_ 21 ^ 
0.5 at these redshifts. 

In order to correctly estimate the contribution of the Lya 
forest to over the entire range of column densities, we 
must adopt a specific model for the dynamics of these clouds. 
Values in excess of the cosmological baryon density of the 
universe derived from nucleosynthesis, Q b hl 0 — 0.05 ± 0.01 
(Walker et al. 1991), would be obtained in the case of spherical 
clouds (Meiksin & Madau 1993; Press & Rybicki 1993). We 
conclude that the large Lya forest clouds must have large 
aspect ratios and contain a large fraction of the baryons in the 
universe. 

2.2. Lyman Limit Systems 

The scarcer Lyman-limit systems (log N H { > 1 7.2), which are 
optically thick to radiation with energy greater than 1 ryd, are 
associated with the halo regions of bright galaxies, and rep- 
resent essentially the same population of absorbers as the 
Mg ii and C iv metal-line systems (Steidel 1993; Stengler- 
Larrea et al. 1995). Steidel (1990) used the column density mea- 
surements of H i and associated ions of heavy elements, in 
conjunction with photoionization models, to estimate the 
physical parameters of these systems. Although there are many 
uncertainties in the analysis, the typical metal abundances are 
Z LLS ~ 0.003 Z Q at redshift z ~ 3, and their mass density 
parameter is 

^lls — 0.008 h 5( } , (3) 

comparable to the mass density of visible matter in nearby 
galaxies. Recent observations obtained with the Keck tele- 
scope have confirmed that the metallicity of a sample of 
Lyman-limit systems with N Hl ~ 3 x 10 1 7 cm 2 is compara- 
ble to that of Lya forest clouds with N Hl > 3 x 10 14 cm 2 
(Cowie et al. 1995). 

2.3. Damped Lyot Absorption Systems 

The damped Lya systems (log N H , > 20.3) may be the pro- 
genitors of present-day luminous galaxies (Wolfe 1990). While 
relatively rare, they account for most of the neutral hydrogen 
seen at high redshifts, with a mass density parameter evolving 
for 0.008 < z < 3.5 as 

n D (z) ~ (3.8 X 10 ” 4 ) exp (0.83 z)h; 0 l , (4) 

for q 0 — 0.5 (Lanzetta, Wolfe, & Turnshek 1995). Thus, at 
z ^ 3, we have n D % 0.005 h^. At redshift z « 2, the typical 
metallicity, Z Dl is approximately 1/10 of solar (Pettini et al. 
1994). A “closed box” cosmic chemical evolution model which 
matches the observed evolution of D^z) also implies that the 
metallicity of the damped Lya systems should decrease sharply 
at z > 3 (Lanzetta et al. 1995). 

2.4. Galaxies 

That star-forming galaxies might rival quasars as sources of 
ionizing radiation is easy to demonstrate. We assume, for pur- 
poses of estimation, that each population is dominated by 
objects near the break of their luminosity function. (This 
approximation is not rigorously true, since a galaxy luminosity 
function with a steep slope at L L* might result in dwarfs 
dominating the luminosity density.) Within factors of 3, the 
Lyc photon luminosities of typical L* disk galaxies and quasars 
are given by S(L*) % 10 53 s" 1 and S(QSO) % 10 56 s -1 . For 
example, the total Lyc production from the Milky Way was 
recently inferred by COBE observations of [N n] 205 pm emis- 
sion to be 3.5 x 10 53 s 1 (Bennett et al. 1994), while the value 
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quoted for QSOs is appropriate for a total luminosity ~ 10 46 
ergs s" K If a fraction </ esc > of the Lyc from O stars escapes the 
galactic H i layer into the IGM, then the ratio of the Lyc 
production rates of L* galaxies to that of QSOs at redshifts 
z = 2-4 is given by 


4>(L*) 1 

S(L*) 1 

U(QSO)J 

S(QSO) 


</esc> * 10</esc> > 


(5) 


where we have adopted comoving densities </>(L*) =10 2 
Mpc -3 and tf>(Q SO) = lO" 6 Mpc 3 . Clearly, if </ esc > * 0.1, 
the two sources are comparable. From models of the escape of 
Galactic Lyc photons from OB associations through “H it 
chimneys” in the H i disk layer, Dove & Shull (1994) conclude 
that this escape fraction is at least 14%. The actual value could 
be somewhat higher if one accounts for the dynamical effects of 
superbubbles in puncturing the H i layer. We return to this 
issue in § 4. 

One problem with the above estimate is that we may not be 
justified in extrapolating the low-redshift density of bright, 
star-forming galaxies to z > 3. In fact, despite considerable 
observational effort, a significant population of primeval, 
young galaxies at high-z has to date evaded detection. While 
the bulk of the B < 24 galaxies are normal dwarfs at low red- 
shifts (Colless et al. 1993), analyses of number counts, colors, 
and the lack of many Lyman-break candidates suggest that 
more than 90% of the faint, blue galaxies with 25 < B < 27 lie 
at z < 3 (Guhathakurta, Tyson, & Majewski 1990; Steidel & 
Hamilton 1993; Madau 1995). To reconcile these observations 
with the hypothesis that galaxies are major contributors to the 
ionizing background at high-z requires that we invoke dwarf 
or low surface brightness galaxies, undetectable to B = 27 at 
z > 3. Such galaxies are likely to have shallow gravitational 
potential wells, thus allowing easy ejection of heavy elements 
produced during their active star-forming phase. 


3. IONIZING PHOTON PRODUCTION 


3.1. The Ionizing Spectrum of Metal-forming Galaxies 

If the metals detected in the Lyot clouds arise from a popu- 
lation of stars with an IMF extending to massive stars (10-100 
M g ), this star formation will be accompanied by the pro- 
duction of hot gas, heavy elements, and ionizing stellar radi- 
ation. Sutherland & Shull (1995) used the evolutionary tracks 
of Schaller et al. (1992), the recently computed heavy element 
yields for massive stars and Type II supernovae (Woosley & 
Weaver 1995), and a recent recalculation of O-star Lyc photon 
production rates (Vacca, Garmany, & Shull 1995) to derive a 
conversion efficiency of rest mass to ionizing continuum, 


^7 Lyc 




M m c 


: 0.002 


( 6 ) 


Here M m is the mass of heavy elements (Z > 6) produced in 
massive stars and L v is the spectral luminosity (ergs s“ 1 Hz *) 
of ionizing photons with energies above hv 0 = 1 ryd. This con- 
version efficiency did not change by more than 10% when the 
IMF slope (T) of stars on the upper main sequence [</N(> M)j 
dM = AM~ r for 8 M 0 < M < 85 Af 0 ] was varied from 
T = 1.0 to 2.0, where T =1.35 corresponds to the usual Salpe- 
ter function (see also Songaila et al. 1990). The increased metal 
yields from high-mass stars are compensated for by a similar 
increase in the Lyc production. 

Figure 1 a shows the time history of the ionizing spectrum 
from one such starburst. Except for the updated OB and Wolf- 



Fig. I. — Modeled input spectrum (Sutherland & Shull 1995) from a Gauss- 
ian starburst of width 2 a = 4 Myr containing 5000 M 0 of stars between 8 and 
85 M 0 with IMF slope T = 1.6 (see § 3.1). Spectrum shows four bands: soft 
UV (5.0-13.6 eV) and ionizing continua of H I (13.6 eV), He I (24.58 eV), and 
He it (54.4 eV). (a) Time history of ionizing spectral luminosity at limes 2, 4, 6, 
12, and 16 Myr (peak is at 6 Myr). (b) Time-integrated spectrum {E v in eq. [9]) 
of the starburst normalized to the total radiated ionizing energy J *£ v dv, 
integrated over the 42.5 Myr lifetime of an 8 M 0 star. 


Rayet stellar spectra, the resulting composite model is similar 
to the young-galaxy spectra of Bruzual & Chariot (1993). The 
stellar He II ionizing continuum above 4 ryd is dominated by 
Wolf-Rayet stars, for which we use the model spectra of 
Schmutz, Leitherer, & Gruenwald (1993). For the purposes of 
this paper, we shall assume that the starburst duration is suffi- 
ciently short compared to a Hubble time that the relevant 
quantity of interest is the time-integrated spectrum (Fig. lb). 
This should provide a reasonable estimate of the ensemble- 
average of a number of starburst populations. 

3.2. The Ultraviolet Background 

The mean specific intensity J v of the radiation background 
at the observed frequency v 0 , as seen by an observer at redshift 

is 

1 f 00 dl (\ + z ) 3 

J(v Q , z 0 ) = T- dz T 77 T u e ( v > z ) ex P [~ T eff K, z 0% z)] , 
4^ Jz a «z ( I + z ) 

(7) 

where e(v, z) is the proper volume emissivity (in ergs s 1 cm 3 
Hz" 1 ) at frequency v = v 0 (l + z)/(l + z 0 ) 9 dl/dz is the line 
element in a Friedmann cosmology, and exp [ — r eff (v 0 , z 0 , z)] is 
the average transmission of a clumpy universe. We have com- 
puted the effective photoelectric optical depth r eff = 
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— ln«e” r » of a cloudy intergalactic medium following 
Madau (1995). For Poisson-distributed clouds we have 
(Paresce, McKee, & Bowyer 1980) 


(V„,z 0 > z)= I dz' | dN m 
Jzq Jo 


8 2 N 

SN Hl d 




(1 


(8) 


where d 2 N/dN H , dz' is the redshift and column density dis- 
tribution of absorbers along the line of sight, and t is the H i 
continuum optical depth through an individual cloud of 
column density N H] . When z 1, z tff becomes equal to the 
mean optical depth. In the opposite limit, the obscuration is 
picket fence-type, and the effective optical depth becomes equal 
to the mean number of optically thick systems along the line of 
sight. The cosmic Lyc opacity is dominated by absorbers with 
t ~ 1. Note that the He i contribution to the attenuation is 
negligible if QSOs contribute substantially to the ionizing 
background, while He n absorption 2 on the way must be 
included when (1 + z) > (1 + z 0 XA 0 / 228 A). 

The ionizing emissivity associated with metal formation can 
be written as 


e < V > z) = n ab., P C ri.(z)C' 2 Z<y esc > . (9) 

where the first term in parentheses is the normalized, time- 
integrated spectrum of our template star-forming galaxy (see 
Fig. lb), fi abs is the baryonic mass density parameter of metal- 
forming material in QSO absorption-line systems, and Z is the 
metal formation rate. To parameterize in simple form the many 
possible evolutionary models for the history of cosmic metal 
formation, we assume that 


Z = const = 


Z(z a ) 

[*h( Z o) — ? H( Z on )] 


( 10 ) 


where t H is the Hubble time, and z on flags the cosmic epoch at 
which metal formation turns on. Equations (7), (9), and (10) 
yield a background flux proportional to the product 
Q abs Z(z 0 )</ esc >. As discussed in § 2, the gaseous baryons in the 
IGM add up to Q abs in the range 0.01 h 50 l to 0.05 h 50 l , while the 
metallicities Z range between 0.01 and 0.1 times the solar value 
of 0.02. Although the product Q abs Z may be highest in the 
damped Lya systems, the Lyc escape fraction from such 
systems may be lower. Note that, at high redshifts, the ionizing 
background will depend only weakly on the assumed metal- 
licity evolution history. At early epochs, sufficient neutral 
hydrogen is contained in the IGM to significantly attenuate 
the Lyc flux from all but the nearby sources, within Az % 0.2. 
Because of this large optical depth, our estimate of the metaga- 
lactic flux at z 0 ss; 3 will depend only on the “local” metal 
formation rate Z. 

Figure 2 shows the 300-1000 A spectrum of the metagalactic 
radiation field from massive stars, computed at z 0 — 3 for dif- 
ferent values of z on == 3.2, 3.5, and 4.0 and filtered through a 
cloudy IGM. We assume a value of n abs Z(z 0 )</ esc > equal to 
l x 10“ 6 , corresponding, for example, to Q abs = 0.02 
(comparable to the baryon density in the Lya forest clouds), 
Z = 2 x 10“ 4 (about 1% of solar, the inferred metallicity of 


2 An approximate (within 5%) formula for the effective H i continuum 
opacity at the observed wavelength < 912 A is z„, z)~0.25x 3 

(x 046 - x° 46 ) + 9.4x f 1 - 5 (.x 018 — x° 18 ) - 0.7x 3 (x‘ ! 32 - x” 132 ) - 0.023 
(x 1 68 - xj 68 ), where x e = (XJ9\2 A)x 0 , x a = (1 + z„), and x ~ (! + z) (Madau 
1995). The first term on the right-hand side represents the contribution of Lya 
clouds, the others are due to Lyman-limit systems. 



Fig. 2. — Observed spectrum at redshift z 9 — 3 of a UV background domi- 
nated by massive stars, modified by the absorption of intervening Lya clouds 
and Lyman-limit systems. We assume a constant metal-formation rate from z on 
to z 0 and adopt n, b# Z(z tf )</ et4: > = 1 x 10" 6 (see § 3.2). The metagalactic flux is 
plotted un units of (Ji 2 0 ) ergs cm -2 s~‘ Hz" 1 sr *. From lop to bottom, the 
three curves correspond to different epochs, z^, = 3.2, 3.6, and 4.0, for the onset 
of heavy-element formation. The spectrum for = 3.2 shows bumpy struc- 
ture because the stellar input spectrum undergoes less redshift smearing. 


the Lya forest and Lyman-limit absorbers), and </ csc ) = 0.25. 
This product is quite uncertain; we return in § 4 to a discussion 
of values of the individual terms. In the three cases considered, 
the background intensity at the hydrogen Lyman edge, J 912 , is 
found to range from 0.13 to 0.36 x 10 21 /if 0 ergs cm 2 s 1 
Hz” 1 sr” 1 . These values are a substantial fraction of direct 
measurements of the metagalactic ionizing flux from the 
“proximity effect” (Bajtlik, Duncan, & Ostriker 1988; 
Bechtold 1994) and are comparable to the estimated contribu- 
tion from observed QSOs (Meiksin & Madau 1993). 

4. DISCUSSION 

Recent observations of the cosmological mass density of 
metals associated with QSO absorption systems at high-z are 
expected to bear significantly on the history of star formation 
and the chemical and ionization evolution of the universe as a 
whole. We have shown that the Lyc photons which accompany 
the production of metals in the Lya forest clouds could con- 
tribute significantly to the ionizing background at high red- 
shifts. In this scenario, the metals are formed from massive 
stars, probably in dense gaseous regions of young galaxies, and 
are then expelled to large distances consistent with the sizes of 
Lya clouds. As a consequence, it is difficult to ascertain the 
exact location of the Lyc production except in an average 
sense. Similarly, the photoionization rates derived here do not 
distinguish between “local” sources (from a single star- 
forming galaxy) and the integrated effects of many such gal- 
axies. One would naturally expect a detailed model to include 
fluctuations in the spectrum due to varying local star- 
formation rates and “filtering” by surrounding H i gas (the 
parameters t, and </ esc > in eq. [9]). 

The model calculations we just described are based on a 
number of parameters, some of which contain significant 
uncertainties. In order to provide an assessment of the validity 
of the stellar photoionization hypothesis, we believe that it is 
important to summarize the weak spots in the argument. As 
shown in equation (9), the Lyc emissivity is proportional to the 
product n abs Z(zJ</ esc >, multiplied by a spectral template and 
an efficiency, rj Lyc % 0.002, with which Lyc radiation is gener- 
ated during the production of heavy elements. The spectral 
template is a fairly robust quantity, since it is based on the 
composite, time-integrated spectrum of massive stars produced 
in a burst of star formation. Likewise, the efficiency rj Lyc is 
nearly constant, as long as the stellar population extends as a 
power law to massive stars (50-80 M G ). The metallicity of the 
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Lya clouds was based on detection of absorption by C iv, and 
not oxygen. Since carbon is produced primarily by low-mass 
stars, there is always the possibility that star formation at high 
redshift may have occurred with an anomalous IMF lacking 
massive stars. We find this possibility slim, but the detection of 
O vi in the Lya forest clouds would settle the issue, since 
oxygen is unambiguously produced in massive stars (Woosley 
& Weaver 1995). 

Next, we address the possible uncertainties associated with 
our estimate of the starburst contribution to the ionizing back- 
ground. We choose to relate these issues to uncertainties in the 
three key terms in equation (9): the baryon density, L2 abs , the 
metallicity, Z(z 0 ), and the Lyc escape fraction, </ esc >. The first 
term is on a fairly solid basis; the baryon density associated 
with all types of Lya absorbers at z % 3 is at least 0.01 of 
closure density and could approach the nucleosynthesis limit, 
Q b » (0.05 + 0.01)fc 50 2 . The Lya forest potentially contains the 
bulk of these baryons, although the uncertainties due to cloud 
ionization and geometry are significant. 

The second term, defining the cloud metallicity history, is 
somewhat more uncertain. Although the metallicities at a 
given epoch are as reliable as the ionization models used to 
derive them, the time history of metal production is less well- 
defined. Current estimates for damped Lya absorbers (Pettini 
et al. 1994) suggest that the bulk of metals formed between 
redshifts z = 2-3; no such estimates exist for the Lya forest. 
For example, the metal-formation rate near the observer red- 
shift z 0 could be negligible, and most of the element production 
may occur at earlier epochs. In this case, the background field 
at z 0 would be quenched by the large continuum optical depth 
of the universe at those redshifts. More likely, galactic chemical 
history, like star-formation history, is spread widely over time. 
In this case, our approximation of a constant metal-formation 
rate may not be too bad. 

The third term, </ esc >, parameterizing the Lyc escape frac- 
tion, is probably the most uncertain of all. This factor depends 
on a variety of stellar, galactic, and interstellar quantities, with 
few observations to guide us. In a theoretical study of the 
structure of H ii regions in vertically stratified galactic H i disk 
layers punctured by OB associations, Dove & Shull (1994) 
used distributions of H i column density (Dickey & Lockman 
1990) appropriate for the solar vicinity, together with the ion- 
izing photon luminosity function of OB associations 
(Kennicutt, Edgar, & Hodge 1989). They concluded that 14% 
of the Lyc photons escaped through the H i layers into the 
halo in the 2.5 kpc region around the solar circle, and they 
speculated on how this fraction might change in other locales. 
For example, the escape of Lyc from a large star-forming 
complex in a galactic nucleus would likely be higher than that 
from a distribution of smaller OB associations in a galactic 
disk. Similarly, the gas distributions and OB formation sites in 
dwarf or irregular galaxies probably differ from those in spiral 
disks. Finally, there is almost no information characterizing 
the structure and types of star-forming galaxies at high red- 
shift. 


A few observational tests might provide some help in con- 
straining </ esc >. Using spectral observations with the Hopkins 
Ultraviolet Telescope, Leitherer et al. (1995) have recently set 
limits of a few percent on the Lyc escaping fraction in low- 
redshift starburst galaxies. Thus, at first sight, the stellar photo- 
ionization hypothesis may be untenable, at least at low 
redshift. However, even though their target galaxies were 
chosen to be bright in the near-ultraviolet, there is no guar- 
antee that the Lyc is not obscured by gas along a different 
orientation (cone effects). In addition, these constraints need to 
be obtained for dwarf galaxies as well. Searching for Ha or Lya 
emission might also constrain the escaping flux, although the 
quantitative sensitivity of this test is likely to be diminished by 
internal dust absorption (Binette et al. 1993). The interpreta- 
tion of Lya emission and Lyc radiation escaping from galaxies 
is intrinsically difficult, since they rely on modeling of optical 
and near-UV radiation to determine the total ionizing radi- 
ation produced. We conclude, therefore, that </ esc > is an 
unconstrained parameter of our model and of provious models 
(Songaila et al. 1990; Miralda-Escude & Ostriker 1990) of 
stellar contributions to the intergalactic radiation field. 

Despite these caveats, we believe that our estimates provide 
some rationale for the hypothesis that hot stars could rival 
quasars as a source of photoionization of the IGM. If this is 
true, the signature of a stellar spectrum should appear in 
several places, notably the He n Gunn-Peterson effect, whose 
detected strength (Jakobsen et al. 1994) requires a relatively 
soft radiation field at 4 ryd compared to 1 ryd (Madau & 
Meiksin 1994; Giroux, Fardal, & Shull 1995). On the other 
hand, if N ( He ii)//V(H i) in the Lya clouds is measured to be 
less than 100, then stars are unlikely to be major contributors 
to the ionizing background; after filtering by the IGM the 
stellar spectra are too steep. Likewise, the distribution of metal 
ionization states of the Lya forest and Lyman-limit absorbers 
should reflect an ionizing stellar contribution, from either local 
sources or the integrated background. 

Therefore, we propose three crucial tests of the stellar photo- 
ionization scenario: (1) confirming the presence of massive 
stars in the high-z stellar IMF by detecting O vi /U1032, 1038 
absorption in the Lya clouds; and (2) constraining the Lyc 
escape, </ esc ), by far-UV continuum and Lya emission obser- 
vations of low-redshift starburst galaxies and by indirect esti- 
mates for star-forming regions in the Galaxy; and (3) 
measuring the ratio N(He ii)/jV(H i) in the lines of sight to 
high-z quasars. 


We thank Ralph Sutherland for his help in deriving the 
“starburst galaxy” ionizing spectrum, and we acknowledge 
helpful conversations on the subject of this paper with Mark 
Giroux, John Stocke, and David Tytler. Support for this work 
was provided by the Colorado Astrophysical Theory Program 
(NASA grant NAGW-766). J. M. S. thanks the Visitor 
Program at the Space Telescope Science Institute where this 
work was completed. 


REFERENCES 


Bajtlik, S., Duncan, R. C., & Ostriker, J. P. 1988, ApJ, 327, 570 
Bechtold, J. 1994, ApJS,91, I 

Bechtold, J., Crotts, A. P. Duncan, R. C., & Fang, Y. 1994, ApJ, 437, L83 

Bechtold, J., Weymann, R. J., Lin, Z., & Malkan, M. A. 1987, ApJ, 315, 180 

Bennett, C. L., et al. 1994, ApJ, 434, 587 

Bergeron, J., & Stasinska, G. 1986, A&A, 169, 1 

Binette, L., Wang, J. C. L., Luo, Z , & Magris, G. 1993, AJ, 105, 797 


Bruzual, G., & Chariot, S. 1993, ApJ, 405, 538 

Colless, M., Ellis, R. S., Broadhurst, T. J., Taylor, K., & Peterson, B. A 1993, 
MNRAS, 261, 19 

Cowie, L. L., Songaila, A., Kim, T.-S., & Hu, E. M. 1995, AJ, 109, 1522 
Dickey, J. M., & Lockman, F. J. 1990, ARA&A, 28, 21 5 
Dinshaw, N., Impey, C. D., Foltz, C. B., Weymann, R. J., & Chaffee, F. H. 1994, 
ApJ, 437, L87 



556 


MADAU & SHULL 


Dove, J. B., & Shull, J. M. 1994, ApJ, 430, 222 
Fan, X.-M. 1994, Ph D. thesis, Columbia Univ. 

Giroux, M., Fardal, M.,& Shull, J. M. 1995, ApJ, 451, 477 
Giroux, M., Sutherland, R. S., & Shull, J. M. 1994, ApJ, 435, L97 
Gruenwald, R., & Viegas, S. M. 1993, ApJ, 415, 534 
Guhathakurta, P., Tyson, J. A., & Majewski, S. R. 1990, ApJ, 357, L9 
Jakobsen, P., Boksenberg, A., Deharveng, J. M., Greenfield, P., Jedrzejewski, 
R., & Paresce, F. 1994, Nature, 370, 35 
Kennicutt, R. C., Edgar, B. K., & Hodge, P. W. 1989, ApJ, 337, 761 
Lanzetla, K. M., Wolfe, A. M., & Turnshek, D. A. 1995, ApJ, 440, 435 
Leitherer, C., Ferguson, H. C., Heckman, T. M., & Lowenthal, J. D. 1995, 
preprint 

Lu, L. M. 1991, ApJ, ApJ, 379, 99 
Madau, P. 1991, ApJ, 376, L33 

. 1992, ApJ, 389, LI 

.1995, ApJ, 441, 18 

Madau, P., & Meiksin, A. 1994, ApJ, 433, L53 
Meiksin, A., & Madau, P. 1993, ApJ, 412, 34 
Miralda-Escude, J., & Ostriker, J. P. 1990, ApJ, 350, 1 

Murdoch, H. S., Hunstead, R. W., Pettini, M., & Blades, J. C. 1986, ApJ, 309, 
19 

Paresce, F., McKee, C., & Bowyer, S. 1980, ApJ, 240, 387 
Petitjean, P., Rauch, M., & Carswell, R. F. 1994, A&A, 291, 29 
Pettini, M., Smith, L. J., Hunstead, R. W., & King, D. L. 1994, ApJ, 426, 79 
Press, W. H., & Rybicki, G. B. 1993, ApJ, 418, 585 
Reimers, D., & Vogel, S. 1 993, A&A, 276, L 1 3 

Sargent, W. L. W., Young, P. J., Boksenberg, A., & Tytler, D. 1980, ApJS, 42, 
41 


Schaller, G., Schaerer, D., Meynet, G., & Maeder, A. 1992, A&AS, 96, 269 
Schmutz, W., Leitherer, C, & Gruenwald, R. 1993, PASP, 104, 1164 
Sciama, D. W. 1993, Modern Cosmology and the Dark Matter Problem 
(Cambridge: Cambridge Univ. Press) 

Songaila, A., Cowie, L. L., & Lilly, S. J. 1990, ApJ, 348, 371 
Steidel, C. C. 1990, ApJS, 74, 37 

. 1993, in The Environment and Evolution of Galaxies, ed. J. M. Shull & 

H. A. Thronson (Dordrecht: Kluwer), 263 
Steidel, C. C., & Hamilton, D. 1993, AJ, 105, 2017 
Steidel, C. C, & Sargent, W. L. W. 1989, ApJ, 343, L33 
Stengler-Larrea, E. A., et al. 1995, ApJ, 444, 64 
Sutherland, R. S., & Shull, J. M. 1995, in preparation 
Tytler, D., & Fan, X.-M. 1994, ApJ, 424, L87 

Tytler, D., Fan, X.-M., Buries, S., Cottrell, L., David, C., Kirkman, D., & Zuo, 
L. 1995, in QSO Absorption Lines, Proc. ESO Workshop, ed. J. Bergeron, 
G. Meylan, & J. Wampler (Heidelberg: Springer), in press 
Vacca, W. D., Garmany, C. D., & Shull, J. M. 1995, ApJ, in press 
Viegas, S. M., & Friaca, A. C. S. 1995, MNRAS, 272, L35 
Viegas, S. M., & Gruenwald, R. 1991, ApJ, 377, 39 

Walker, T. P., Steigman, G., Schramm, D. N., Olive, K. A., & Kang, H. 1991, 
ApJ, 376, 51 

Weymann, R. J. 1993, in The Environment and Evolution of Galaxies, ed. J. M. 

Shull, & H. A. Thronson, (Dordrecht: Kluwer), 213 
Wolfe, A. M. 1990, in The Interstellar Medium in Galaxies, ed. H. A. Thronson 
& J. M. Shull (Dordrecht : Kluwer), 387 
Woosley, S., & Weaver, T. 1995, ApJS, 101, 181 



